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We report here the results of a study of a partial sample of 2.773 x 10 7 K — > 7r 7r°7r° decays 
recorded in 2003, showing an anomaly in the 7r°7r° invariant mass (Moo) distribution in the 
region around Moo = 2m+ where m+ is the charged pion mass. This anomaly has never 
been seen in earlier experiments. It can be interpreted as an effect of the charge exchanging 
scattering process ty*ty~ — > 7r°7r°. 



1 Beams, detectors and reconstruction 

During the years 2003-04, the experiment NA48 at CERN SPS has collected a sample of fully 
reconstructed ~ 2 x 10 8 K — ► 7r ± 7r 7r decays in order to study direct CP violation by compar- 
ing the Dalitz plot distributions of K + and K~ . The experiment makes use of two overlapping 
simultaneous focused kaon beams of opposite charge and momentum of 60 GeV/c ±3.8%, se- 
lected by a system of "achromat" magnets and collimators. The two beams enter in a 114 m 
long decay volume and the tt^tt -^ final state is reconstructed by combining the signals com- 
ing from a spectrometer and from a liquid Krypton calorimeter (LKr)^. The spectrometer, 
consisting of four drift chambers^ and a dipole magnet located between the second and the 
third chamber, allows to track the charged pion and to measure its momentum. The resolution 
is cr(p)/p = 1.02% © 0.044%p, where p is in GeV. The calorimeter is used to reconstruct the 
vr° -> 77 decays. The resolution on 7 energy is a{E)/E = 0.032/Vi? © 0.09/£ © 0.0042. The 
space resolution on the transverse coordinates can be parametrized as a x = a y = 0.42/\/^©0.06 
cm (E in GeV). Events with at least one charged particle track and at least four energy clusters 
in the LKr are selected for further analysis. The distance between any two 7-rays in the LKr 
is required to be larger than 10 cm and, in addition, the distance between each 7-ray and the 
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Figure 1: (Left) Invariant mass distribution of Tr*^ ^ candidate events. (Middle) 7r°7r° invariant mass squared 
distribution. (Right) Zoom of the previous in the region around the value Moo = 2m + indicated by the arrow. 

impact point of any track on LKr must exceed 15cm. Other cuts ensure full containment of the 
electromagnetic shower in the LKr. The constraint of a common decay vertex for the two tt° 
is used to pair the four photons. From now on the correct pairing is assumed to be (1,2) and 
(3,4). Momentum conservation and small angle approximation allow to compute the distance D 
of the decay vertex to the LKr and the two 7r° invariant mass Moo as: 

D = 2k (^i E 2 rh + jEkE trj^ 
Moo = (l/D)xjY*~^Erf j W 

where Ei, Ej are the energies of the i-th and j-th 7-ray and rij is their transverse distance at the 
LKr plane, mo is the 7r° mass. Note that both D and Moo computation involve only quantities 
measured by the LKr and no charge track parameters are used. This fact implies that even 
in the case of ir^ — ► ^ decays in flight, Moo is correctly computed. Figure Jflleft) shows the 
invariant mass distribution of the system ir^it^Tr® . The distribution is dominated by the 
peak, as expected. Tails originate from tt^ — > decays in flight. The selection requires that 
the vr^vr ^ invariant mass differs from the nominal mass quoted in the PDG by no more 
than ±6 MeV. The fraction of events with wrong 7-ray pairing is ~ 0.25% as estimated by a 
Montecarlo simulation. 

2 7r°7r° invariant mass distribution 

Figure (^middle) and (fright) show the two neutral pions invariant squared mass Mqo- The 
r.m.s of the Mq resolution curve increases with Mq , varying between ~ 0.0002 GeV 2 and 
~ 0.001 GeV 2 at the end of the allowed Mq range. This excellent resolution is the result of 
the intrinsic energy and spatial resolution of LKr. A sudden change in the slope (Cusp) near 
Moo = 2m + can be clearly seen. At this point the resolution on Moo is ~ 0.0003 GeV 2 . Such 
an anomaly has not been observed in previous experiments. 

3 Interpretation of the cusp 

The observed sudden change of the slope suggests the presence of a threshold "Cusp" effect 
from the decay — ► tt^^tt^ contributing to K — ► tt^^-k amplitude through the charge 
exchange reaction tt + tt~ — > 7r°7r°. This phenomenon has been recently discussed by Cabibbo^ 
who computed the — > tt^tt ^ amplitude taking into account the 1-loop diagram shown in 
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Figure 2: 1 loop diagram contributing to K — > h ± -k°'k . 

figure [2 The diagram produces a discontinuity when the looping charged pions pair flips from 
off to on mass shell at the value Moo = 2m + . The PDG parametrization for the — > tt^'k^-k 
decay amplitude M.q is given by 



Mo = ^ (I + -guf + h'u 2 (2) 

u = (S3 — SQ)/m^_ is the Lorentz-invariant variable where Sj = (Pk — Pi) 2 (i = 1,2,3), Pk 
and Pi are the 4-momentum vectors respectively of the initial Kaon and of the three outgoing 
pions; so = (s\ + S2 + sz)/3 and of course Mq = S3. In the Cabibbo theory, the diagram in 
figure [21 is responsible for a new term in the amplitude Ai% oc (do — (12) J proportional to the 
difference between the 1 = and 1 = 2 S-wave irir scattering lengths. The term J changes from 
real to imaginary at Moo = (2m+). The distructive interference of Mq and .Mi in the total 
amplitude M. = M.Q + M\ is responsible for the cusp and for the apparent lack of event below 
the threshold. More recently Cabibbo and Isidori™ have extended the amplitude calculation at 
2-loops level suggesting a precise way to extract the (do — 02) parameter from the data. The 
most recent theoretical prediction (do — a2)m+ = 0.265 ± 0.004 is given in reference^. 

4 Fit and results 



Model X 2 /ndf 

a la PDG (full region) 15000/139 

a la PDG (above cusp) 240/198 

Cabibbo 1-loop 463/149 

Cabibbo-Isidori 2-loops 159/147 

2-loops and pionium 154/146 



Table 1: fit against various models 



Table^shows the x 2 obtained when the experimental Mq spectrum is fitted against various 
model. The naive PDG description cannot account for the cusp structure observed and it gives 
a "crazy" % 2 although it seems to fit reasonably above the cusp. The Cabibbo 1-loop theory is 
still not adequate to fit the spectrum while the 2-loop calculation is quite satisfactory. 

FigureElat left shows the best fit to data in the interval 0.074 < Moo < 0.097 GeV 2 obtained 
by slightly modifying the Cabibbo-Isidori model to account for a small excess of event located 
exactly at Moo = 2m+ on top of the observed cusp. This excess can be interpreted as formation 
of pionium. The pionium is the electromagnetic bound state 7r + 7r~ which is expected to decay 
mainly in 7r°7r° with a time constant 10 -16 sec. Figure 01 at right shows the distribution of the 
difference between data and fit when the pionium contribution is set to zero in the fit. A fixed 
amount of pionium derived by the estimation published in® has been used for the final fit. The 
value found for the scattering length parameter is (do — d2)?n + = 0.281 ±0.007 where the error is 
only statistical. A very conservative estimation of the systematic uncertainties is given in table 




Figure 3: (Left) Best fit to data. (Right) Excess of event at Moo = 2m+ interpreted as pionium formation. 

Source Error 
skipping pionium ±0.008 
Cut on photon to track distance ±0.004 
Dependence on vertex position ±0.009 

K+/K- difference ±0.006 

total ±0.014 

Table 2: Systematics error on (do — a2)m+ 

12 We tried to exclude 10 points around the cusp from the fit to reduce the sensitivity to the 
pionium component. Other checks have been performed by changing the cut on the minimum 
distance between photons and charged pion at the LKR or by selecting events in different region 
of the decay vertex. As a preliminary result we quote the value 

(oo - a 2 )m + = 0.281 ± 0.007(stat.) ± 0.0U(syst.) ± 0.0U(theo.) 

The theoretical error, equal to the 5% of the central value, has been suggested by Cabibbo and 
Isidori as the achieved accuracy in their model. 

5 Conclusions 

The 7r°7r° invariant mass (Moo) measured from a sample of 2.773 x 10 7 K — > tt^tt ^ fully 
reconstructed decays collected by NA48/2 experiment at the CERN-SPS shows an anomaly at 
Moo = 2m+. This anomaly can be explained by a simple model based on contribution from the 
decay — > tt ± tt + 'k~ through the charge exchange reaction tt + it~ — > 7r°7r°. A high-statistics 
measurement of the Moo distribution could be used therefore to provide a precise determination 
of the 7T7T scattering parameter (ao — 02)- The expected increase of the event sample statistics 
by more than a factor 5 coming from the data collected in 2004 should provide a substantial 
reduction of the conservative systematics quoted here. The quality of the data calls for an 
improved theory including isospin breaking effects and radiative corrections. 
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